home *** CD-ROM | disk | FTP | other *** search
/ IRIX 6.2 Development Libraries / SGI IRIX 6.2 Development Libraries.iso / dist / complib.idb / usr / share / catman / p_man / cat3 / complib / COMLR2.z / COMLR2
Text File  |  1996-03-14  |  5KB  |  133 lines
  1.  
  2.  
  3.  
  4. ____CCCCOOOOMMMMLLLLRRRR2222((((3333FFFF))))                                                        ____CCCCOOOOMMMMLLLLRRRR2222((((3333FFFF))))
  5.  
  6.  
  7.  
  8. NNNNAAAAMMMMEEEE
  9.      COMLR2, SCOMLR2  -  EISPACK routine.  This subroutine finds the
  10.      eigenvalues and eigenvectors of a COMPLEX UPPER Hessenberg matrix by the
  11.      modified LR method.  The eigenvectors of a COMPLEX GENERAL matrix can
  12.      also be found if  COMHES  has been used to reduce this general matrix to
  13.      Hessenberg form.
  14.  
  15.  
  16. SSSSYYYYNNNNOOOOPPPPSSSSYYYYSSSS
  17.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee  ccccoooommmmllllrrrr2222((((nnnnmmmm,,,, nnnn,,,, lllloooowwww,,,, iiiigggghhhh,,,, iiiinnnntttt,,,, hhhhrrrr,,,, hhhhiiii,,,, wwwwrrrr,,,, wwwwiiii,,,, zzzzrrrr,,,, zzzziiii,,,,
  18.      iiiieeeerrrrrrrr))))
  19.           iiiinnnntttteeeeggggeeeerrrr           nnnnmmmm,,,, nnnn,,,, lllloooowwww,,,, iiiigggghhhh,,,, iiiieeeerrrrrrrr
  20.           iiiinnnntttteeeeggggeeeerrrr           iiiinnnntttt((((iiiigggghhhh))))
  21.           ddddoooouuuubbbblllleeee pppprrrreeeecccciiiissssiiiioooonnnn  hhhhrrrr((((nnnnmmmm,,,,nnnn)))),,,,hhhhiiii((((nnnnmmmm,,,,nnnn)))),,,,wwwwrrrr((((nnnn)))),,,,wwwwiiii((((nnnn)))),,,,zzzzrrrr((((nnnnmmmm,,,,nnnn)))),,,,zzzziiii((((nnnnmmmm,,,,nnnn))))
  22.  
  23.           ssssuuuubbbbrrrroooouuuuttttiiiinnnneeee ssssccccoooommmmllllrrrr2222((((nnnnmmmm,,,, nnnn,,,, lllloooowwww,,,, iiiigggghhhh,,,, iiiinnnntttt,,,, hhhhrrrr,,,, hhhhiiii,,,, wwwwrrrr,,,, wwwwiiii,,,, zzzzrrrr,,,, zzzziiii,,,,
  24.      iiiieeeerrrrrrrr))))
  25.           iiiinnnntttteeeeggggeeeerrrr           nnnnmmmm,,,, nnnn,,,, lllloooowwww,,,, iiiigggghhhh,,,, iiiieeeerrrrrrrr
  26.           iiiinnnntttteeeeggggeeeerrrr           iiiinnnntttt((((iiiigggghhhh))))
  27.           rrrreeeeaaaallll              hhhhrrrr((((nnnnmmmm,,,,nnnn)))),,,,hhhhiiii((((nnnnmmmm,,,,nnnn)))),,,,wwwwrrrr((((nnnn)))),,,,wwwwiiii((((nnnn)))),,,,zzzzrrrr((((nnnnmmmm,,,,nnnn)))),,,,zzzziiii((((nnnnmmmm,,,,nnnn))))
  28.  
  29.  
  30.  
  31. DDDDEEEESSSSCCCCRRRRIIIIPPPPTTTTIIIIOOOONNNN
  32.      On INPUT
  33.  
  34.      NNNNMMMM must be set to the row dimension of two-dimensional array parameters
  35.      as declared in the calling program dimension statement.
  36.  
  37.      NNNN is the order of the matrix.
  38.  
  39.      LLLLOOOOWWWW and IGH are integers determined by the balancing subroutine  CBAL.
  40.      If  CBAL  has not been used, set LOW=1, IGH=N.
  41.  
  42.      IIIINNNNTTTT contains information on the rows and columns interchanged in the
  43.      reduction by  COMHES, if performed.  Only elements LOW through IGH are
  44.      used.  If the eigenvectors of the HESSEN- BERG matrix are desired, set
  45.      INT(J)=J for these elements.
  46.  
  47.      HHHHRRRR and HI contain the real and imaginary parts, respectively, of the
  48.      complex upper Hessenberg matrix.  Their lower triangles below the
  49.      subdiagonal contain the multipliers which were used in the reduction by
  50.      COMHES, if performed.  If the eigenvectors of the Hessenberg matrix are
  51.      desired, these elements must be set to zero.  On OUTPUT The upper
  52.      Hessenberg portions of HR and HI have been destroyed, but the location
  53.      HR(1,1) contains the norm of the triangularized matrix.
  54.  
  55.      WWWWRRRR and WI contain the real and imaginary parts, respectively, of the
  56.      eigenvalues.  If an error exit is made, the eigenvalues should be correct
  57.      for indices IERR+1,...,N.
  58.  
  59.      ZZZZRRRR and ZI contain the real and imaginary parts, respectively, of the
  60.  
  61.  
  62.  
  63.                                                                         PPPPaaaaggggeeee 1111
  64.  
  65.  
  66.  
  67.  
  68.  
  69.  
  70. ____CCCCOOOOMMMMLLLLRRRR2222((((3333FFFF))))                                                        ____CCCCOOOOMMMMLLLLRRRR2222((((3333FFFF))))
  71.  
  72.  
  73.  
  74.      eigenvectors.  The eigenvectors are unnormalized.  If an error exit is
  75.      made, none of the eigenvectors has been found.
  76.  
  77.      IIIIEEEERRRRRRRR is set to Zero       for normal return, J          if the J-th
  78.      eigenvalue has not been
  79.         determined after a total of 30*N iterations.  Calls CSROOT for complex
  80.      square root.  Calls CDIV for complex division.  Questions and comments
  81.      should be directed to B. S. Garbow, APPLIED MATHEMATICS DIVISION, ARGONNE
  82.      NATIONAL LABORATORY
  83.  
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  
  92.  
  93.  
  94.  
  95.  
  96.  
  97.  
  98.  
  99.  
  100.  
  101.  
  102.  
  103.  
  104.  
  105.  
  106.  
  107.  
  108.  
  109.  
  110.  
  111.  
  112.  
  113.  
  114.  
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  
  121.  
  122.  
  123.  
  124.  
  125.  
  126.  
  127.  
  128.  
  129.                                                                         PPPPaaaaggggeeee 2222
  130.  
  131.  
  132.  
  133.